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ABSTRACT 

We show that the small scatter around the Fundamental Plane (FP) of 



(N 
> 

oo : massive elliptical galaxies can be used to derive important properties of their 

\D . 

£NJ i dark and luminous matter. The central velocity dispersion cro, appearing in 

OO ■ 

' ( e -g-) the Fundamental Plane, is linked to photometric, dynamical and geo- 



metrical properties of (luminous and dark) matter. We find that, inside the 
effective radius R ei the matter traced by the light must largely dominate over 
the dark matter (DM) , in order to keep the ellipticals close enough to the FP. 



2 ' This recalls analogous findings for spiral galaxies 



6 
b 

In particular we also find that cuspy DM distributions, as predicted by 

numerical simulations in ACDM cosmology, are unable to explain the very 
• i— i . 

. existence of the FP; in fact, according to this theory, the structural proper- 

■ ties of dark and luminous matter are so interwoven that a curved surface is 

predicted in the log-space (cto, R e , L), rather than a plane. In order to agree 
with the FP, CDM halos must have concentrations parameters in the range 
of 5 — 9 (i.e. values significantly lower than the current predictions). 

Assuming a more heuristic approach and allowing for cored DM halos, we 
find that the small intrinsic scatter of the FP yields to i) an average value 
for the dark-to-light-traced mass ratio inside the length-scale of light R e of 
about 0.3, ii) a mass-to-light ratio of the matter traced by the light increasing 
with spheroid luminosity: M sp h/L r oc L® 2 in Gunn-r band, with a value of 



5.3 at = 2.T x 10^ 
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2 Borriello, Salucci & Danese 
1 INTRODUCTION 

In the hierarchical scenario, dark matter (DM) halos have driven, from a variety of ini- 
tial conditions, a dissipative infall of baryons and formed the galactic systems we observe 
today (White and Rees, 1978). Thus, we expect that DM halos exist within and surround- 
ing any galaxy, regardless of its luminosity and morphological type. This prediction had 
overwhelming confirms for disk galaxies, due to the existence of good dynamical tracers 
and their intrinsic simple geometry (see Persic and Salucci, 1997). Elliptical galaxies (E's), 
however, are much more complicated objects, due to their 3-dimensional shape, stellar or- 
bital structure and velocity dispersion anisotropy. These factors have made ambiguous the 
interpretation of observational data. 

A number of different mass tracers have been used to probe the gravitational potential in 
tenth of E's and derive their mass distribution: integrated stellar absorption spectra, X-ray 
emission from hot gas, rotating gas disks, motions of globular clusters or satellite galaxies 
and, in last years, weak gravitational lensing. As result, the presence of dark matter in 
E's, especially in the external regions ( £ 10 kpc), is proven (e.g. Loewenstein and White, 
1999). On the other hand, a kinematical modeling of the inner regions (i.e. within the half- 
luminosity "effective" radius R e ), has been performed for only a small number of ellipticals 
(e.g. van der Marel, 1991; Saglia et al. 1992, 1993; Bertin et al, 1994; Kronawitter et al, 2000; 
Gerhard et al., 2001); the results point to a tendency for moderate dark matter amounts 
inside R e . 

Since its discovery, the "Fundamental Plane" (Djorgovski & Davis, 1987; Dressier et 
al., 1987) has been one of the main tools to investigate E's properties: effective radius R e , 
central velocity dispersion <jq and mean effective surface brightness J e of spheroidal galaxies 
are linearly related in the logarithmic space and galaxies closely cluster on a plane, with a 
surprisingly low orthogonal scatter. To explain these linear relations between photometric 
and dynamical quantities in log-space, most studies on the Fundamental Plane (FP) have 
considered models in which the mass is distributed parallel to light. However, in presence 
of non-baryonic dark matter, this hypothesis is an obvious oversimplification and, at least, 
unjustified. Indeed, this would a priori require either: i) dark and luminous component are 
distributed according the same profile, thus revealing a similarity of properties and behavior 
which seems very unlikely or ii) the dark matter component is always negligible with respect 
to the luminous matter. 
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Dark Matter in Elliptical Galaxies 3 
Within the above framework, in this paper we address the following issues: 

• to derive the relation between the central velocity dispersion o"o and the mass distri- 
bution parameters, including the effect of a dark matter halo. In particular, we assume a 
spherical model with an isotropic luminous component and a dark halo, more diffuse than 
the spheroid, 

• to reproduce the observed Fundamental Plane and, therefore, to constrain the mass 
distribution in E's, 

• to discuss the results in the light of Cold Dark Matter predictions. 

Considering elliptical galaxies as two-components systems, complementary strategies are 
possible. One chooses a distribution function for both components and then imposes specific 
constraints from the observations. The other includes the ordinary stellar component (or, 
better, any traced by light (TBL) mass component) in a frozen spherical halo. The former 
approach is helpful in exploring the self-consistency of the dynamical configuration (e.g. 
Ciotti, 1999). The latter, we will adopt in this paper, has the advantage of providing a 
simpler connection between observational quantities and the parameters of the mass model. 

The outline of this paper is the following: in §2 we describe two-components models, 
whose mass distributions are shown in §3. In §4 we derive and discuss the velocity dispersion 
(line-of-sight profile and central value) predicted by the mass models we consider. In §5, we 
introduce the data, fit the models to the Fundamental Plane and discuss the results. Finally, 
conclusions are presented in §6. Throughout the following work, we assume, where needed, 
a flat ACDM Universe, with Q m = 0.3, Q A = 0.7, h = 0.7 and <j 8 = 1.0. 

2 THE VELOCITY DISPERSION IN THE 2-COMPONENTS MASS 
MODELS 

The observed velocity dispersion and, in particular, the central velocity dispersion <jq is a 
fundamental dynamical property related to the gravitational potential of both the dark halo 
and the TBL component (sometime, we call the "luminous stellar spheroid" by this way, 
to recall the possible existence of a number of BH's and/or an amount of non-baryonic 
DM, perfectly mixed with the luminous stellar matter). Let us start by assuming a spherical 
and non-rotating stellar system, with the stellar velocity dispersion which is the same in 
all directions perpendicular to a given radial vector. If af(r) denotes the stellar velocity 
dispersion along the radial vector and cr^(r) the dispersion in the perpendicular directions, 
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4 Borriello, Salucci & Danese 

the Jeans hydro dynamic equation for the mass density traced by the light p sp h{f) in the 
radial direction reads (Binney and Tremaine, 1987): 
dp sph (r)a 2 (r) 2 (3(r)p sph (r)a 2 (r) _ GM(r ) 

dr + r r 2 Psph[r) [ } 

with the boundary condition p sp h(r)a 2 (r) — > for r — > oo. In eq.([T]), the parameter /3(r) 

1 — cr|(r) / a 2 (r) describes the anisotropy degree of the velocity dispersion at each point, with 
P = 1,0, — oo for completely radial, isotropic and circular orbit distributions, respectively. 

Dynamical analysis of ellipticals exclude substantial amount of tangential anisotropy 
and find (3( < R e ) ~ 0.1 - 0.2 (e.g. Matthias and Gerhard, 1999; Gerhard et al, 2001; 
Koopmans and Treu, 2002) for objects of different luminosities. For reasons of simplicity, 
then, in calculating the central velocity dispersion to be used for statistical studies over a 
large sample of galaxies, we assume an isotropic velocity dispersion tensor: (3 = 0. 

Eq.(P connects the spatial velocity dispersion of the component traced by the light to 
its density profile and to the total matter distribution M(r) = M sp h(r) + M^(r). Under the 
hypothesis of isotropy, the above equation assumes the well-known integral form, in which 
we single out the halo term: 

2/ \ g r Ps P h{r')M{r') , _ 2 2 

^r W = / — 72 dr = a ns P h( r ) + °r-A r ) (2) 

Psph{r) Jr r' 2 

As external observers of galaxies, we measure only projected quantities. Let R be the 
projected radius and the surface stellar mass density. As usual, we take into account 

that mass density and the spatial velocity dispersion are related to the surface mass density 
E(i2) and to the projected velocity dispersion o~p(R) by the two Abel integral equations for 
the quantity p sp h and p sp hO~ 2 . Then, a second step, consisting in a further integration along 
the line of sight, allows us to obtain the (observed) velocity dispersion profile o~p(r): 

a2 ^ R) = m E p *X>-* r dr s ahph{R) + aUR) (3) 

where E{R) = Jr [2 r p sph {r)/{r 2 - R 2 ) 1 ' 2 } dr. 

As spectro-photometric observations are performed through an aperture, let us define 

o~a{Ra) as the luminosity-weighted average of op within a circular aperture of radius Ra- 
2tt r R A 

°a(Ra) = I o- 2 P {R) I(R) R dR = a% sph (R A ) + a% h (R A ) (4) 

where I(R) is the surface brightness profile I{R) = S(i2)/T (assuming the stellar mass- 
to-light ratio T constant with radius) and L(Ra) = 2 tt J^ a I(R) R dR is the aperture 
luminosity. 
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Dark Matter in Elliptical Galaxies 5 

The dynamical quantity in the Fundamental Plane is the "central" velocity dispersion o~ , 
which, observationally, corresponds to the projected velocity dispersion luminosity-weigthed 
within the aperture of the observations. The velocity dispersion data are to be brought to a 
common system, independent of the telescope and galaxy distance: this is done by correcting 
them to the same aperture of R e /8 (J0rgensen et al., 1996), which is typical of measurements 
of nearby galaxies. Therefore, we compare model and observations by calculating o~o as the 
luminosity-weighted o~p{R) within R& = 1/8 R e . 

The resulting velocity dispersion profiles, o>(r), crp(R) and cta(Ra), can be all expressed 
as the sum of two terms (eqs. 0, [3] and f|): the first one is due to the self-gravity of the 
spheroid traced by the star light (labelled by sph); the second one (labelled by h) is due to 
the effect of the luminous-dark matter gravitational interaction and, therefore, it charges 
relevance according the characteristics of the DM distribution. 



3 THE MASS DISTRIBUTION 

3.1 The distribution of the mass traced by light 

We describe the component traced by the star light by means of a Hernquist (1990) spher- 
ical density distribution, that is a good approximation to the de Vaucouleurs R 1 ^ law (de 
Vaucouleurs, 1948) when projected and, at the same time, allows analytical calculations: 

/ \ M S pfo k R e 

PsAr) = 17r(r + U) s (5) 
where M sp h is the total mass traced by the star light and k ~ 0.55. The mass profile derived 

from eq.@ is: 

M.Ar) = M, rkJ J0 w (6) 

The Hernquist functional form, of course, cannot reproduce the fine features of the surface 
brightness profile (e.g. boxy isophotes, small variations in slope); neverthless, it is sufficient 
for our aims, since we will just consider large scale properties in the mass distribution of 
objects belonging to a large E's sample. 

3.2 The DM distribution: ACDM halos 

N-body simulations of hierarchical collapse and merging of CDM halos have shown that 
gravity, starting from scale-free initial conditions, produces an universal density profile that, 
for r — > 0, varies with radius as r~ a , with a ~ 1 — 2 (e.g. Navarro, Frenk & White, 1997, 
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6 Borriello, Salucci & Danese 

hereafter NFW; Fukushige and Makino, 1997; Moore et al, 1998; Ghigna et al, 2001), 
weakly dependent on the cosmological model. We adopt the well-known NFW halo density 
profile: 

pNFw{r) = . , — (7) 
(r/r s )(l + r/r s ) 2 

where r s is the inner characteristic length-scale, corresponding to the radius where the 
logarithmic slope of the profile is —2. It results convenient to write the NFW mass profile 
as: 

where A(x, y) = ln(l +x/y) — (1 +y/x)~ 1 for any pair of variables (x, y). The concentration 
parameter is defined as c = r vir /r s ; r vir and M vir are, respectively, the halo virial radius 
and mass. The definition of the virial radius is strictly within the framework of the standard 
dissipationless spherical collapse model (SCM); however, also in more realistic hierarchical 
models, it provides a measure of the boundary of virialized region of halos (Cole and Lacey, 
1996). 

Considering a dark halo at redshift z, the virialized region is the sphere within which 
the mean density is A vir (z) times the background universal density at that redshift (pbkg = 
p c (1 + z) 3 , with p c the critical density for closure at z — 0). The virial mass is defined as: 
M vir = | 7r A V i r (z) pbkg rl ir , with the virial overdensity A V i r being a function both of the 
cosmological model and the redshift: for the family of flat cosmologies (fi m + ^A = 1 ) , it can 
be approximated by (Bryan and Norman, 1998): A vir (z) ~ 18vr 2 + 82(Q(z) - 1) - 39(Q(z) - 
l) 2 /il(z). From the above equations, we derive: 

r vir = 0.142 A mr (z)~ 1/3 (1 + z)- 1 ^ kpc (9) 

A fundamental result of CDM theory is that the halo concentration well correlates with 
the virial mass: low-mass halos are denser and more concentrated than high-mass halos 
(Bullock et al., 2001; Cheng and Wu, 2001; Wechsler et al., 2002) in that, in average, they 
collapsed when the Universe was denser. Numerical experiments by Wechsler et al. (2002) 
for a population of halos identified at z = show that: 

with en ~ 20.8. The Poisson error for galactic halos in the mass range M vir ^ 1O 11_13 M is 
less than 10% and, virtually, no halo is found with c < 12. 
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Dark Matter in Elliptical Galaxies 7 

This is the first mass model we consider here: it is composed by a stellar bulge with a 
Hernquist profile embedded in a spherical dark NFW halo (hereafter H+NFW model). We 
neglect the effects of a possible adiabatic coupling between baryons and dark matter because 
it has the effect of increasing the DM density inside R e by adiabatic compression and, then, 
of worsening the fit to the data (see § |5.3.1[ ). The "total" dark-to-TBL mass ratio, defined as 



T V i r = M V i r /M sp h is a crucial parameter of the mass model. Both BBN predictions about the 
primordial DM-to-baryons ratio and CMB anisotropy observations point at a lower limit 



for T vir of ~ 



3.3 The DM distribution: cored halos 



In last years, studies of high resolution rotation curves of spiral and dwarf galaxies casted 
doubts on the presence of the central cusps predicted by cosmological simulations of DM 
halos. Actually, the observational results suggest that dark halos are more diffuse than the 
luminous component and their densities flatten at small radii, whereas the stellar distribution 
peaks towards the centre (e.g. Moore, 1994; Flores and Primack, 1994; Burkert, 1995; de 
Battista and Sellwood, 1998; Salucci and Burkert, 2000; de Block, McGaugh and Rubin, 
2001; Bordello and Salucci, 2001). 

An useful analytic form for halos with soft cores has been proposed by Burkert (1995) 
for dwarf galaxies and, then, was extended to the whole family of spirals by Salucci and 
Burkert (2000): 

pB{r) = (1 + r/r )[l + (r/r ) 2 } (H) 
This profile is characterized by a density-core of extension r and value po, while it resembles 

the NFW profile at large radii. From eq.(|ll|) the mass profile reads: 

w^-H (12) 

where, for any pair of variables (x,y), B(x,y) = — 2 arctan(x/7/)+2 ln(l+x/y)+ln[l + (x/y) 2 } 
and M e is the dark mass within R e . Then, in analogy of spiral galaxies, we propose a mass 
model consisting of a Hernquist bulge plus a Burkert dark halo (hereafter H+B model). 
This is the second mass model we investigate in this paper; let us notice that, in this case, 
the halo mass distribution is characterized by two free parameters: the dark-to-TBL mass 
ratio within the effective radius T e = Mdm I 'M ap h \r b and the halo core radius in units of the 
effective radius r /R e . We assume r £ R e to ensure a constant DM density in the region 
where the stars reside: otherwise the model would essentially coincide with the NFW one. 
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Figure 1. Line-of-sight velocity dispersion due to the spheroid (in units of total velocity dispersion) as a function of radius 
for: left) H+NFW case (different values of R e and T v i r are indicated); right) H+B case (different values of ro/R e and T e are 
indicated). 

It is worth noticing that the Burkert profile is purely heuristic and featureless in the 
optical regions of galaxies; as a consequence, unless data at large radii are available, it is not 
possible to determine the halo virial radius and the total galaxy mass. 

4 MASS-VELOCITY DISPERSIONS RELATIONS FOR CUSPY AND 
CORED MODELS 

We compute the velocity dispersion profiles, including the effect of a spherical dark halo, 
for both H+NFW and H+B cases. We resolve eq.(Q), eq.(|3|) and eq.(|]) by assuming the 
density/mass profiles of eq.(|5]), eq.(|G]) and eq.(|^) in H+NFW case and of eq.(|5]), eq.(|G]) and 
eq.(|T2"D in the H+B one. The detailed calculus in given in Appendix. 

In Fig.|l| (left) we show, for H+NFW models, the radial profile of the quantity crp h /ap, 
the line-of-sight velocity dispersion due to the TBL component, in units of the total l.o.s. 
velocity dispersion. We take M sp h = 2 x 1O 11 M ; however, the mass dependence is very weak 
and the curves in Fig.|l| (left) are well representative of those with stellar masses in the range 
~ 5 x 10 9 — 1O 12 M . We consider different plausible values for the total dark-to-TBL mass 
ratio T vir and R e . Let us notice that, once we fix the total halo mass M vir = T vir ■ M sp h, the 
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Figure 2. The functions i^NFW an <i entering the expression of ao, for H+NFW model (left) and H+B model (right). 

halo characteristic radius r s is completely determined via the c—M vir relationship: therefore, 
different curves in Fig.|I] (left) correspond to 2 & r s /R e & 30. 

We realize that the contribution of a CDM halo to the velocity dispersion can be large 
(~ 50%) even at small radii R ^ R e /2, in galaxies with large effective radii and/or small 
values of r s /R e , indipendently of the value of T vir . 

In Fig.|l] (right) we plot the TBL-to-total l.o.s. velocity dispersion ratio for H+B models, 
for different values of the parameters r /R e and r e . Notice that the profiles just depend on 
these parameters and it is not necessary to assume specific values both for M sp h or R e . The 
main consequence of the smooth halo profile in the H+B model is that the halo contribution 
to <Jp(R) is low at small R, even for models with a relevant amount of DM in the central 
region (e.g. (T 2 P . sph /a 2 P ~ 80% at -R e /3 when Y e = 1). The velocity dispersion in the central 
regions is more directly connected to the properties of the TBL mass distribution. 

The relationship of the spheroid mass with central velocity dispersion tT and the effective 
radius R e depends on the DM mass distribution. Recalling that er = aA(R e /8), we find: 



a 2 = (0.174 + r„ ir F NFW ) 



G M. 



sph 



R f 



a 2 = (0.174 + T e F B ] 



G M 



sph 



Re 



H + NFW 



H + B 



(13) 



(14) 



where -F NFW and F B are shown in Fig.^[ F NF w depends on r s (M vir )/ R e and, very weakly, 
on M vir = T vir ■ M sp h. The term in the r.h.s. of eqs.(13) and (14) shows that, according to 
the values of the model parameters, <jq could be strongly affected by the DM gravitational 
potential, so that the derivation of the spheroidal mass from obsevational properties critically 
depends on the actual DM mass profile. This is shown in (Fig.^, left), where we plot the 
TBL-to-total erg, assuming M sp h = 2 x 1O 11 M and different model parameters: the TBL 
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Figure 3. The TBL mass component contribution to ag for different models parameters: left) H+NFW predictions for M sp h = 
2 X 10 11 Mg, as function of r„; r and for different r s /R e ; right) H+B predictions, as function of the parameters T e and ro/Re- 

contribution to a 2 , is under dominant also at small r s /R e and becomes almost negligible for 

30. 

In H+B mass models, Fb weakly depends on ro/R e , so that we can assume Fb — 3.6 x 
10 -2 ; the TBL contribution to ctq (Fig.^|, right) remains the dominant one, even for an 
amount of DM within R e comparable to the luminous one. Moreover, from eq.(|14|) we infer 
that, for cored configurations (i.e. with r > R e ), a is weakly dependent on the DM internal 
amount (Fb is of the order of 10~ 2 ): as a matter of fact, it just increases of ~ 30% when 
T e varies of a factor 3. This is a natural consequence of the smoothness of the dark matter 
distribution with respect to the more concentrated distribution of the luminous spheroid. 



5 FITTING MASS MODELS TO THE FUNDAMENTAL PLANE 
5.1 The Sample 

We build the data sample from several works by J0rgensen, Franx & Kjaergaard (hereafter 
JFK). They provide spectroscopy and multicolour CCD surface photometry of E/SO galaxies 
in nearby clusters. The photometric data are from JFK (1992) and (1995a) in Gunn-r, 
their passband with the largest quantity of data. The spectroscopic measurements are taken 
from JFK (1995b) and references therein. Out of the whole JFK sample, we selected a 
homogeneous subsample of 221 E/SO galaxies in 9 clusters, including Coma, whose properties 
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Figure 4. Histograms of sample galaxies properties: (a) effective radius R e , (b) cllipticity e at R e , (c) central velocity dispersion 
(To an d (d) luminosity L r in Gunn— r band . 

are shown in electronic form in Tab.l at the URL: www.sissa.it/ap/ftp. In particular, we 
rejected spiral, interacting, peculiar and field galaxies (due to the greater uncertainty of 
their distance). For each cluster, we adopt the distance derived in JFK (1996). The FP r.m.s. 
scatter is 0.084 in logi? e (this is equivalent to a ~ 17% uncertainty in galaxy distances that, 
then, could be the main scatter contributor). Typical measurement errors are Alogi? e = 
±0.045, AlogJ e = ±0.064, Alogcr = ±0.036, and AlogL r = ±0.036, maybe large enough 
to imply that the FP is free from an intrinsic scatter. 

The statistical distributions of effective radius, ellipticity at R e , observed central velocity 
dispersion and Gunn-r luminosity of the selected galaxies are shown in Fig.^. It is worth 
noticing that the sample galaxies are distributed around = 2.7 x lO lo L r0 , the character- 
istic luminosity of the ellipticals luminosity function in r-band (Blanton et al., 2001) and 
that most of the objects have little/moderate ellipticity (< e >= 0.29 ± 0.17) and then, a 
reasonably spherical stellar distribution. 



5.2 Forcing models to the Fundamental Plane 

The physical interpretation of the Fundamental Plane assumes the virial theorem to be the 
main constraint to the structure of ellipticals. Assuming elliptical galaxies to be 1-component 
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log R e 




Figure 5. The surface in the log-space (Be, L r , o"o) where we expect galaxies in the H+NFW scenario for Y v i r = 30, compared 
to sample data (red points). The units are: R e in kpc, L r in L r © and ctq in km/s. 

and homologous systems, the virial theorem and the existence of the FP imply a slow but 
systematic variation of the mass-to-light ratio M/L with the luminosity, whose physical 
origin is debated. This homology also determines a quasi-linearity of the relations connecting 
the gravitational and kinematic scale parameters of the galaxies to the observables o"o and 
R e (see Prugniel and Simien, 1997). However, according to the properties of the second, 
dark, mass component, this property could be lost, in such a way that the gravitational and 
photometric scales are not anymore connected in a simple, log-linear way. 

We adjust the mass models parameters to fit the observations in the log-coordinates 
space: effective radius R e , central velocity dispersion cr = aA(R e /8) an d total luminosity 
in Gunn-r band, defined as L = 2nR 2 e I e . The effective surface brightness I e in L /pc 2 is 
calculated from \i e in mag arsec -2 : log/ e = — 0.4/x e — 26.4, for Gunn-r band (JFK, 1995a). 
In fitting the log-surface cr (R e , L r ) to the observations, we leave free the mass-to-light 
ratio M sp h/L r of the TBL component and, respectively, F v i r in the H+NFW case and T e in 
the H+B case. In the latter, we assign a constant value to the parameter r§jR e = 2, similar 
to results for spirals (Borriello and Salucci, 2001), as the fit depends very weakly on it. We 
characterize the mass-to-light ratio as T r = Msph/Lr = T r * (L r /L*) a , with Y r * and a free 
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Figure 6. The mass model parameters, obtained by best-fitting the H+NFW surface in the log-space (R e , L r , <rrj) to the FP 
(68%, 95% and 99% CL). 

parameters. Here, we neglect a possible weak dependence of T r on R e , but we will discuss 
this point later. 

5.3 Results and discussion 

5. 3. 1 H+NFW mass model 

This model is unable to provide a plane surface in the log-space (00, R e , L r ), for plausible 
values of the free parameters. In Fig.[5| we show the effect on the surface by adding to the 
stellar spheroid a dark NFW component with the reasonable dark-to-TBL mass ratio of 
T«r = 30. The surface curvature prevents us from properly fitting the data, especially in 
the region occupied by galaxies with large effective radius and low luminosity, for which the 
DM contribution to o~ is unacceptably high. 

The results of the fitting procedure are shown in Fig.|6|, with contours representing 68%, 
95% and 99% CL. The best-fit model is consistent with no DM at all (T V i r = 2 ± 4 at 1 a) 
and values of Y vir £ 10 are excluded at > 95% CL. Solutions very marginally permitted 
(T vir ~ 10 — 15) still require a very high efficiency of collapse of baryons in stars (~ 90%). 
This is at strong variance with the inferred budget of the cosmic baryons (Salucci and Persic, 
1999) and with current ideas of galaxy formation, for which feedback mechanisms (such as 
SN explosions and central QSO activity) transfer thermal energy into the ISM and inhibit 
an efficient star formation (e.g. Dekel and Silk, 1986; Romano et al., 2002). They predict 
T vir to be significantly higher than the cosmological value of ~ 8. 
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Figure 7. Concentration parameter en and T r * (68%, 95% and 99% CL) in agreement with the observed FP, for H+NFW 
mass model with V V i r = 30. 

The mass-to-light ratio of the TBL component in Gunn-r band is found ~ 5, i.e. well 
within values predicted by passive evolution of an old stellar population, generated with a 
standard IMF (e.g. Trager et al., 2000). As a matter of fact, this encourages us to consider 
the TBL mass as composed of just stellar populations. 

The curvature of the surface in the logarithmic space of the observables is a conse- 
quence of the particular concentration-mass relation predicted by CDM N-body simula- 
tions. Anyway, once we assume a value for T vir , we can consider different concentration 
parameters (for example, by moving a fraction of DM outside R e ). In order to illustrate 
the consequences, we have explored the realistic case of T vir = 30. For simplicity, we as- 
sume c = c\\ ■ (M„j r ,/1O 11 M )~ 013 , and we look for values of en for which CDM profiles 
are in agreement with the Fundamental Plane. In Fig.[7] we show the value for c\\ implied 
by the narrowness of the observed FP. In order to have a small DM fraction inside R e and 
recover the FP, we must lower the concentration parameter down to ~ 5, well below the 
standard predictions of numerical simulations of halos in ACDM cosmology (e.g. Wechsler 
et al., 2002). Again, it is worth noticing the robustness of the estimate of the mass-to-light 
ratio of the TBL component, which results to be insensitive to even such a model change. 

5.3.2 H+B mass model 

In this case, the presence of dark matter, distributed independently of the stellar profile, 
does not alter the FP surface shape, which still remains a plane. The best-fit mass model is 
obtained for (see Fig.|): 
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Figure 8. The H+B plane (edge-on), best fitting the data. 
£ x 0.21±0.03 



(15) 
(16) 



(at 68% CL). In Fig|| (left) we show the 68%, 95% and 99% confidence contours for the 
TBL mass-to-light ratio parameters Y r * and a. In the fit, the parameter T e is somewhat 
correlated to T r *: the greater DM amount in the bulge region, the lower the mass-to-light 
ratio of the TBL mass component. In Fig.|9| (right) we show this correlation, marking the 
T e range corresponding to 99% CL in Y r *. Notice that a variation of T e of a factor ~ 5 
corresponds to a much smaller variation of T r *, giving prominence to the strong stability of 
the stellar mass-to-light ratio we deduce from fits. 

Finally, to check the reliability of the fit against our assumption of spherical stellar 
distribution, we perform the model fit by considering only 133 galaxies with small ellipticity 
(e < 0.3); the resulting best-fit parameters are consistent with those of the whole sample, 
with a difference in the mean values of about ~ 5%. 

It is worth stressing that assuming T e constant in the fit is not the most general possibility 
and, in principle, the DM contribution within R e can well vary with luminosity. Although 
the investigation of this point is beyond the possibilities of our database, we notice that 
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Figure 9. Best-fit H+B mass model: left) 68%, 95% and 99% CL for the stellar mass-to-light parameters; right) the correlation 
between T r * and r e , the dark-to— stellar mass ratio within R e . Dashed lines mark the 99% C.L. in T r » and T e . 

a variation would be in agreement even with recent dynamical and photometric studies of 
mass distribution in individual ellipticals (Kronawitter et al., 2000; Gerhard et al., 2001). 
Moreover, Gerhard et al. (2001) found M DM /M tot \ Re ~ 10 — 40%; in comparison we find 
M DM /M tot \ Re ~ 30%. 

From eq.flilj) we derive the relation between a\R e jG and the spheroid mass or, equiva- 
lently, the total mass within the effective radius M e ~ (1 + r e ) 0.42 M sp h: 



M, 



sph 



G 



M e = M sph (R e ) + M h (R e ) ~ 2.9 



erg R e 
G 



(17) 
(18) 

0.29. Eq.flT|) explicitly 



using Fb = 3.6 x 10~ 2 (see Fig. 2) and the best-fit value T e 
contains the effect of a DM halo and is to be compared with the "gravitational" mass at 
R e (M G ~ 2 o"q R e /G), found by Burstein et al. (1997), taking the standard Keplerian 
formula M e = R e V^ ot /G and assuming Vj? ot = 3<7g. Eg. (|17D , instead, is in good agreement 
with results by Ciotti, Lanzoni and Renzini (1996); indeed they find for their HP model (a 



mass configuration similar to our H+B model) M sp h = cm 



with cm — 3 — 6, according 



to the value of the total dark-to-TBL mass ratio (in the range 10 — 70). 



In Fig.ITU] we show the distribution of the mass-to-light ratio in Gunn-r band of the 
TBL component, obtained by inserting in eq.(|l7]) the observed <To, R e and L r . Continuous 
line is the mean correlation provided by the FP fit : 



Tf r — 5.3 (yL r I 



\0.21 



(19) 
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Figure 10. The distribution of the stellar mass-to— light ratio in Gunn-r band. The continous line is the mean correlation 
provided by the FP fit. The typical data error at 1 a is also shown. 

By testing the residuals of the mass-to-light ratio as function of the effective radius R e , we 
find no correlation within the statistical errors: T r oc 00±0 - 05 . A possible weak dependence 
on the effective radius, therefore, seems not sufficient to justify the scatter observed in the 
luminosity dependence of M sph /L. 

Since part of the galaxy sample has also been observed in different photometric bands 
(JFK, 1992; JFK, 1995a; see Tab.l), we investigated the M sp h/L variations with luminos- 
ity (for a smaller number of galaxies) in Johnson U and B and Gunn-f band, obtaining, 
respectively, the slopes 0.34 ± 0.05, 0.27 ± 0.03 and 0.25 ± 0.06, with similar large scatter, 
thus independent of the photometric band. The slope of the relation M sp h/LB oc L° B 27±om 
seems to be smaller than the value of ~ 0.6 ± 0.1, obtained from velocity dispersion profiles 
analysis (Gerhard et al., 2001). Anyway, this slope tends to higher values when we con- 
sider a DM fraction (within R e ) decreasing with galaxy luminosity. For example, we obtain 
M sph /L B oc L^ 45±0 05 assuming T e = 1 for L B < 10 10 L QB and V e = 0.1 for L B > 10 10 L QB . 

Recent observations have enlightened the potential of galaxy-galaxy (weak) lensing to 
probe the DM halo around galaxies at large radii, where it has been impossible, so far, to 
find kinematic tracers. Studies of weak lensing by SDSS collaborations (McKay et al., 2002; 
Guzik and Seljak, 2002) find M 2 eo/L r ~ 110, for a L rJf elliptical galaxy with a NFW dark 
matter halo (M 2 qo is the mass projected within an aperture of radius 260 h^ 1 kpc). Since 
the first data bin is at R = 75 ft, -1 kpc, the SDSS g-g lensing is not sensitive to small scales 
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Figure 11. Properties of H+B mass model for a L* galaxy (T r = 5.3, T e = 0.3): left) total (dark+TBL) mass-to— light ratio 
in Gunn— r band vs. radius, for different values of halo core radius; right) dark-to— TBL mass ratio as function of the radius, 
for different values of ro/R e . 

(where NFW and Burkert profiles actually differ); therefore, we can compare our results 
with SDSS ones. We estimate, at large radii, the dark-to-TBL mass ratio and the (total) 
mass-to-light ratio of the H+B mass model. In Fig.|TT] (left) we show the total (dark+TBL) 
mass-to-light ratio for a L r * elliptical galaxy. Assuming R e ~ 5 — 10 kpc for L = L r *, 
the 260 h^ 1 kpc aperture corresponds to ~ (35 — 70)R e , for h = 0.7. We obtain the value 
M/L r ~ 110 at this aperture for halo core radii in the range ~ (2.4 — 2.8) R e . This value 
for the core is in agreement with the findings in spirals (Borriello and salucci, 2001) and 
suggests some form of DM-baryons interplay as the origin of the soft density cores in halos. 

In Fig.|TT] (right) we show the cumulative dark-to-luminous mass ratio M^/Msph, the 
value we found above for the halo core radius of tq ~ (2.4 — 2.8)R e corresponds to a mass 
ratio, at very large radius, of Mh/M sp h — 15 — 30. Remarkably, this range of values is in good 
agreement with results from spectro-photometric models of E's formation in the spheroid 
mass range 3 x 10 9 ^ M sph ^2x 10 11 (Romano et al., 2002; their Fig.10), models including 
even chimical evolution and feedback. 
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6 SUMMARY AND CONCLUSIONS 

We have shown that the very low scatter of elliptical galaxies around the Fundamental 
Plane can be statistically used to put very interesting constraints on DM distribution within 
them. The central velocity dispersion is the key quantity we have dealt with. We have briefly 
reviewed its relationship with the mass distribution of both traced-by-light and dark matter. 
Then, we selected a sample of 221 E/SO galaxies with L r > 2 x 10 9 L Q in 9 clusters, endowed 
with very good photometric and spectroscopic data. The sample defines the classical FP in 
the log-space (a , R e , L r ), with the expected small scatter (0.084 in log R e , to be compared 
to a measurement uncertainty Alogi? e = ±0.045). 

We tested the reference model of cuspy DM distribution, namely the NFW model, and the 
cored model proposed by Burkert (1995). Our analysis shows that these luminous galaxies 
are largely dominated within the effective radius by matter traced by light, independently 
of the DM distribution model, cuspy or cored. In particular, for the cuspy NFW model, we 
have shown that the small scatter of our sample galaxies around the Fundamental Plane 
severely challenges the ACDM predictions. In such a theory, the structural properties of 
dark and luminous matter are so interwoven that in the log-space (o"o, R e , L) they produce 
a curved surface, rather than a plane, for plausible values of the total dark-to-TBL mass 
ratio. We conclude that in order to keep the small scatter around the FP we have either 
to keep F vir unacceptably low or to decrease the halo concentration well below the value 
currently predicted by simulations in ACDM cosmology. 

Considering a cored DM density distribution, the agreement with the observed FP im- 
plies a dark-to-luminous mass fraction within the effective radius of ~ 30% and a lu- 
minosity dependence of the spheroid mass-to-light ratio in Gunn-r band: M sp h/L r = 
(5.3 ± 0.1)(L r /L ;(tr ) ' 21±0 ' 03 . An important result is the robustness of the mass-to-light ratio 
of the spheroidal component we obtained, which is in good agreement with predictions by 
stellar evolution models. 

It is also worth noticing that, besides for spiral and dwarf galaxies, a cored DM halo, 
with low internal (within 2 — 3R e ) density which increases as r~ 3 at larger radii, is successful 
to explain also the structure of elliptical galaxies, pointing to an intriguing homogeneous 
scenario. Within this framework, we argue that dark matter in E's can be investigated by 
a reasonably large number of galaxies with measures of l.o.s. velocity dispersion at ~ R e . 
Although, so far, such observations have been severely hampered by the steep decreasing of 
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the surface brigthness with radius, higher and higher sensitivity reached by recent surveys 
offers a good view to obtain a better resolution of the two mass components, in the whole 
region where baryons reside. 
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APPENDIX A: VELOCITY DISPERSION IN DETAIL 

In this appendix, we will detail the procedure to compute velocity dispersions by means of 
Jeans hydrodynamic equations. Let us set all radii in units of the effective radius: f = r/R e , 
R = R/R e , r = r /R e , f s = r s /R e and f vir = r vir /R e . 

Mass distributions 

Traced-by-light mass component - the radial profiles of mass density, mass and surface mass 
density are (Hernquist, 1990): 

PsAr) = ^ Flif) (A1) 

M sph (r) = M sph F 2 (f) (A2) 

E <*> = 2^ ^ (A3> 
where k ~ 0.5509 and: 

F '< f > = TJfTW (A4) 

™ = WTW (A5) 

F >w = ] < A6 > 

with X(R) = [1 - (R/k) 2 ]- l l 2 §ec\r\R/k) for < R < k 
and 

X(R) = [{R/kf - lJ-^Sec- 1 (R/k) for R > k. 
Dark matter halo - the NFW mass profile reads: 

M NFW {r) = M sph F vir (A7) 

where, for any pair of variables (x,y), A(x,y) = ln(l + x/y) — x/{x + y). In particular, 
recalling that c = r vir /r s , we have A(r vir , f s ) = ln(l + c) — c/(l + c). 
The Burkert halo mass profile is: 

M B (f) = 0.416 M sph T e ||^| (A8) 
where B(x, y) = — arctan(x/|/) + 2 ln(l + x/y) + ln[l + (x/y) 2 }. 
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By inserting eqs. Q), (Q, Q£3|), (|£7|) and O in eqs. <g), (|) and (§ (after variables 
substitutions), we will obtain the velocity dispersions profiles: 
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Spheroid self— interaction terms 

Out of the stellar spheroid self-interaction terms in velocity dispersions, o~ 2 . sph , o~ 2 P . sph and 
c\. h , the first two can be analitically obtained (Hernquist, 1990): 



GM, 



a r;sph\ r ) 



sph 



12 k R P 



12 . , ff + k" 

— r(r + k) In I 



r + k 



52 42 „ 12 . 
25 + — r + — r 2 + — r 
k k k 



a 



where: 



J_ GM sph F 4 (R) 
6 k 



R e 

RT\ 
k 2 J 



-3 



F 3 (R) 
_ R 2 



(A9) 
(A10) 



R G 



k 2 X(i?) 8 k* 



p4 p2 \ 

28 ^ + 35 p - 20 + 



i? 6 R 4 R 2 ' 
- 24 F + 68 ^~ 65 P +6 



O 7T — 

k 



(All) 



We obtain the luminosity weigthed velocity dispersion in the aperture Ra = Ra/Rc, by 
integrating eq.( |A10| ): 

1 GM sph J Ra F 4 (R) R dR 



a A;sph( R A) 



(A12) 



6 k R e jRa F3 (jfyR d R 

where integrals must be numerically performed. For the aperture R a — 1/8, we obtain the 
stellar contribution to the "central" velocity dispersion: 
1 GM sph Jo /8 F 4 (R) R dR 



2 

a 0;sph 



6 k R e fV 8 F 3 (R)RdR 
GM sph 



0.174 



Re 



(A13) 



Spheroid— halo interaction terms 

We apply the same procedure for calculating the luminous-dark matter interaction terms: 
o~ 2 . h , a 2 P . h and o\. h . In this case, however, the integrations are always numerical. 



H+NFW mass model 
'Msph _ 

R e A(f vir ,f s ) Fi(f 



a 2 NFW (r) = A ,--, k T^h — (A14) 



poo Fl (r') A(f',f s ) , , 

Jb r — — = dr 



al;NFw(k) = 2 _ ^ (A15 ) 
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df 



df 



GM h r Io Ra R dR 

a\. NFW (R A ) = 2 k 3 Z—ZfEt j 1 (A16) 

' Re A(r vir ,r s ) J RA F 3 (R)RdR 

where we can reduce the free parameters to the only virial mass i) by using eq.([|) for the 

virial radius, with A vir (z = 0) ~ 337: 

rw ~ 2.59 x 10- 2 f^Y" kpc (A17) 

and ii) by assuming ACDM c — M„ r correlation (Wechsler et al., 2002, their Fig. 16), which, 
together with eq. (|A17| ), gives: 



{M mr ) fi / M, 



0.48 



. mry ^ VUJ kpc (Ai8) 

c(M mr ) V 10 M Q y K J 



Burkert halo: 



a r 2 . B (r~) = 0.416 — ^— £ — (A19) 



J'OO ~ Jf. 'pi ^~ 

a 2 (£) = 0.831 k 3 — ^— (A20) 

Re 5(1, f ) ^(A) 



TOO Fl (f') B(f',f ) , 

<B (^) = 0.831 *■ ^ "~ Sg= jg (A21) 

5(1, r ) J Ra F 3 (R) R dR 

From eqs. ( |A16j ) and ( |A21| ), by setting i? a = 1/8, we obtain the halos contributions to the 



central velocity dispersions given in eq.(^) and fll4}) , where the functions F NFW and F B are 
defined as: 

k 3 1 



J 1/8 F 3 (R) RalR A(f vir (M vir ),f s (M vir )) 



roc F^f') A(f' ,f a (M„ ir ) , 

1/8 JR r ~T dr 

^ ^ 

o F 3 (R) 1 ; 

F t (f') B(r',f ) , , 
TOO ~ Jf JF2 *™ 7~ 

-Fb ^ 0.83 -r^ . . . — . . / R di? (A23) 

Jo /8 F 3 (R) RalR B(l,r ) Jo F 3 (R) 

with the constant of proportionality: k 3 / J^ 8 F 3 (R) R dR ~ 8.31. Let us notice that F^fw 
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depends on both R e and M vir = T vir ■ M sph , while F B is only function of the parameter 
f = r /R e . 



© 0000 RAS, MNRAS 000, 000-000 



